Quantum Hall Ferromagnetism in Graphene 



00 

o 
o 

(N 



Kentaro Nomura^ and Allan H. MacDonald^ 

^Department of Physics, University of Texas at Austin, Austin TX 78712-1081, USA 

(Dated: May 12, 2008) 

Graphene is a two-dimensional carbon material with a honeycomb lattice and Dirac-like low- 
energy excitations. When Zeeman and spin-orbit interactions are neglected its Landau levels are 
four- fold degenerate, explaining the 4e^//i separation between quantized Hall conductivity values 
seen in recent experiments. In this paper we derive a criterion for the occurrence of interaction- 
driven quantum Hall effects near intermediate integer values of e^//i due to charge gaps in broken 
symmetry states. 
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Introduction — Two-dimensional graphite (graphene) is 
a gapless semiconductor with a honeycomb lattice and 
an unusual massless Dirac-fermion band structure 1^ that 
has long attracted theoretical attention. The topology of 
its Bloch states leads to large momentum-space Berry 
phases[2!|, quantized and half-quantized Hall effects, and 
a vanishing density of states at the neutral Fermi en- 
ergy which qualitatively alters the way in which electron- 
electron interactions [3] influence electronic properties. 
The integer quantum Hall effect in graphene is expected 
to be unusual because its Landau levels are widely sepa- 
rated and fourfold degenerate in the absence of weak Zee- 
man and spin-orbit interactions. Interest in graphene has 
increased recently because of experimental progress y, 
including the discovery of the integer quantum Hall 
effect [5| with quantized values of the Hall conductivity 
{axy = 4(n + l/2)(e^/ft-)) separated by 4e^//i. In this 
Letter we address the quantum Hall effects that should 
occur at intermediate integer values of filling factor v, 
giving rise to plateaus at axy = v{e'^/h), in principle 
for all integer values of v. These additional plateaus are 
expected to arise from charge gaps induced by electron- 
electron interactions, but have not yet been observed. 
They would be a new example [6| of the enhanced inter- 
action physics that occurs at integer filling factors in a 
strong magnetic field whenever N > 2 Landau levels are 
degenerate or nearly degenerate. At integer filling fac- 
tors the mean-field-theory scenario in which symmetries 
are broken to open gaps between quasiparticle orbitals 
usually applies. The ground state is then well approxi- 
mated by an unrestricted Hartree-Fock state [7| in which 
an integer number i < N oi Landau levels associated 
with orthogonal SU{N) spinors is occupied. The phe- 
nomenon of interaction induced gaps and broken symme- 
tries at integer filling factors is known as quantum Hall 
ferromagnetism. The four-fold degeneracy of graphene's 
Landau levels follows from approximate spin-degeneracy 
and from Bloch state degeneracy between two inequiva- 
lent points in the honeycomb lattice Brillouin zone. The 
low-energy physics of graphene is well described [3| in a 
four-component spinor envelope-function formalism with 
SU{A) invariant electron-electron interactions. Graphene 



is therefore a good example of SU{A) quantum Hall 
ferromagnetism (8[, much more accurately approximating 
this symmetry than bilayer electron systems fyl for exam- 
ple. The absence of additional integer Hall plateaus due 
to quantum Hall ferromagnetism in all but the most re- 
cent samples is almost certainly due to disorder. Fig. [T] 
summarizes the estimate of the minimum sample mo- 
bility required to see quantum Hall ferromagnetism in 
graphene which is explained below. 

Massless Dirac- Weyl quasiparticles — The k ■ p Hamilto- 
nian of the graphene bands is 



Hq = v{px + eAx)Tzcrx + v{py + eAy)(jy 



(1) 



where Tz — ± labels the two-degenerate [K and K') val- 
leys, (Ta are Pauli matrices that act in the space of the 
two-atoms per unit cell, and A{r) is the vector poten- 
tial. In the zero field case, the Hamiltonian ^ has lin- 
ear dispersion E = ±vhk for both spin states and for 
both K and K' valleys. In a magnetic field the spec- 
trum of Hq consists of four-fold degenerate (including 
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FIG. 1: Phase Diagram for 5(7(4) quantum Hall ferromag- 
netism in the n = and n = 1 Landau levels of graphene. 
In our model the ordered region is bounded by a maximum 
value of Vs, the ratio of the density of Coulomb scatterers 
to the density of a full Landau level. Vs is inversely propor- 
tional to the product of the sample mobility and the external 
field strength and order near integer filling factors requires the 
minimum values for this product indicated on the right-hand 
vertical axis. 



spin) Landau level branches with £"„ = ±.hvy/2n/ i as in- 
dicated in the inset in Fig.l. For n = eigenfunctions 
of different valleys are localized on different honeycomb 
sublattices while for n f^ 0, they are symmetric or anti- 
symmetric combinations of two-dimensional-electron-gas 
Landau level n states on one sublattice and level n + 1 
states on the other sublattice. The spin degeneracy is 
lifted by weak Zeeman coupling that we neglect for the 
moment. 

Stoner Criterion — The simplest approximation [10] for 
interacting electrons in a strong magnetic field is one 
in which interactions are treated in the Hartree-Fock 
approximation and disorder in the self-consistent Born 
approximation,lQ,](SCBA). In the strong field limit the 
total energy in this approximation is given by 
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where N^ is the orbital Landau level degeneracy, fia- is 
the Fermi level for spinor-component cr, 



X = 



(27r)2 



V.ntiq) eM-Q^^V^) F\q) (3) 



is the exchange integral, £ = {hc/eBY^'^ is the magnetic 
length, and F{q) is a form factor we will discuss later. 
In the SCBA the Landau level spectral function has the 
form 



A{E) = — [l-{E/Tff'\ 
where F is the Landau level width|10l|: 
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|C/rf.,(g)|2 exp(-g2^V2)F2(q) (5) 



which is estimated below. In Eqs.Q and ^ Vint{<l) 
and Udisio) are the Fourier transforms of the electron- 
electron interaction and the disorder potential and rig 
is the density of disorder scatterers. In the case of 
graphene the spinor index a runs over four possible val- 
ues. These expressions assume perfect S'[/(4) invariance 
of the disorder-scattering and electron-electron interac- 
tions. While this is certainly an approximation, we be- 
lieve it to be an excellent one. They also assume that 
spatial invariance is recovered after disorder averaging, 
so that the electron density matrix is diagonal in its or- 
bital labels and the energy simply proportional to the 
number of orbitals in a Landau level Ncf,. In Eq.([2|) the 
i>a values are the eigenvalues of the density matrix in 
spinor space. 



dE A{E), 



(6) 



which are invariant under unitary transformations of 
the four-dimensional spinor-space, consistent with SU{4) 
symmetry. 



In the normal state the four (nearly) degenerate Lan- 
dau levels are equally occupied. To look for broken sym- 
metry states we write v^ ^ ^ + Sv^ where v^ — N/N^ is 
the total filling factor in the four-fold degenerate Landau 
level of interest. Expanding to second order in Sva- and 
using that ^^ Si^a = wc find that 
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where /io is the normal state Fermi level. The normal 
state is unstable when the second term in square brack- 
ets is negative, in other words when XA{^q) > 1. This 
criterion for quantum Hall ferromagnetism is closely anal- 
ogous to the Stoner lll| criterion for ferromagnetism in 
metals, and has been successfully applied [12] to under- 
stand the appearance of spin-splittings at odd integer 
filling factors in a semiconductor two-dimensional elec- 
tron gas. In the case of quantum Hall ferromagnetism 
(QHF) the exchange energy competes with disorder en- 
ergy rather than with band energy. We can apply the 
QHF Stoner criterion to graphene by relating the dis- 
order potential to the zero-field mobility of graphene, a 
quantity that is conveniently available from experiment. 
Zero-Field Mobility and Coulomb Scattering — We start 
from the Boltzmann transport theory expression for the 
conductivity, applied to the four-fold degenerate Bloch 
bands of graphene: 
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where t is the scattering rate. 
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(9) 

9 is the scattering angle, kp is the Fermi wavevector 
and q — 2kp sin(6'/2) is the scattering wavevector on 
the circular two-dimensional Fermi surface. The last 
6'-dependent factor in Eq.( ^ is non-standard and is 
due to the wavevector dependence of the relative phase 
of graphene Bloch band wavefunctions on the two sites 
within its honeycomb lattice unit cell. The factor of 
kp in Eq.( ^ reflects the density dependence of the 
density-of-final states for elastic scattering of Fermi sur- 
face quasiparticles. For short-range scatterers the inte- 
gral in Eq.([9|) remains finite as density ex kp vanishes. 
Since 'D{Ep) is proportional to kp for two-dimensional 
Dirac bands, Eq.([9l) implies a conductivity that is inde- 
pendent of kp and therefore independent of carrier den- 
sity. Indeed theoretical studies of the conductivity of 
graphene [13[ predict that the conductivity has a weak 
density dependence, remaining finite as kp — > 0. Ex- 
periment, on the other hand, finds that the mobility 



/i = —a/ne in graphene is nearly constant except at very 
low-densities and that it has values ~ 10''cm^V~^s~^ in 
samples that are sufficiently high quality to exhibit the 
integer quantum Hall effect. Evidently quasiparticle scat- 
tering amplitudes are enhanced at lower densities in such 
a way as to convert the kp^ dependence of the scatter- 
ing rate in Eq.([9]) to a k^^ dependence. One plausible 
explanation for this behavior is that Dirac band quasi- 
particle scattering is dominated by Coulomb scattering 
from charged defects near the graphene plane. For two- 
dimensional graphene Udisiq) = Vc{q) = 27re^/g. Insert- 
ing this expression in Eq. ( [9]) we obtain that 
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where g = |- ^ 3 is the effective fine structure constant 
used to characterize the ratio of Coulomb interaction and 
band energy scales in graphene. In Eg.f llOp us should be 
thought of as the density of Coulomb scatterers that are 
located in the substrate within a Fermi wavelength of the 
graphene plane. The influence of more remote scatterers 
is suppressed by the factor exjp(—qd) that appears in the 
two-dimensional Fourier transform of the Coulomb inter- 
action. Inserting Eq.([TO| in Eq.([8|) we find that mobility 
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In systems with Coulomb electron-electron or electron- 
impurity interactions screening normally plays an es- 
sential role, changing long-range interactions into short- 
range ones. In a static approximation, the screened dis- 
order potential in graphene is 
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where Ti{q) is the polarization function of the graphene 
Dirac bands. Screening does not change the density de- 
pendence of the conductivity in graphene because U{q) 
also scales like kp. The influence of screening on the 
mobility can be estimated by making a Tomas-Fermi 
approximation, replacing Il{q) by U{q = 0) = 'D{Ef)- 
When the coupling constant g is much larger than 1 
Usc{q) — {hv'K)/{2kp) and 
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yielding a value for the mobility that is 16g^ times larger 
than the unscreened value. We note that g cancels in 
Eq.([T31), which is fortunate because its effective value 
can be influenced by non-universal substrate dielectric 
screening. Corrections to Eq. ( [13]) becomes impotant for 
g < 1. We can use these expressions to extract a value 
for the density of scatterers n^ from measured mobili- 
ties. This procedure might retain partial validity, de- 
pending on the details, even if the limiting scatterers 



are not Coulombic. Similar density dependence could 
in principle arise from a partially accidental combina- 
tion of disorder sources that gives rise to a similar in- 
crease in transition rates at small wavevectors, or from 
strong short-ranged scattering that approaches the uni- 
tary limit. Other potential disorder sources include ran- 
dom crumpling of the graphene sheet and coordination 
defects in the graphene sheet that give rise to long range 
strain fields. The procedure we now use to translate be- 
tween zero-field mobilities and strong-field Landau level 
widths will retain its validity in some, but not all, plausi- 
ble scenarios. In particular the values of iiB at the stoner 
phase boundary are likely to be similar for Coulomb and 
topological defect scatterng[lj|. 

Self- Consistent Screening in a Magnetic Field — We are 
now in a position to estimate the Landau level width and 
apply the Stoner criterion. For Coulomb scattering 
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Notice that F^ diverges if we neglect screening. We 
now need to specify the form factor F{q). Taking the 
Coulomb interaction to be diagonal in honeycomb lattice 
site index it follows jlOf that the form factor F{q) = 1 for 
n = and that 

Fiq) = I [L|„|(<z2£V2) + L|„|_i(<z2^V2)] (15) 

for n ^ 0. If the magnetic field is strong enough to ne- 
glect coupling between different Landau levels the normal 
state polarization function Il{q)\ is given approximately 
by 
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The factor of 4 in Eq. ( [TB)) is the graphene Landau level 
degeneracy and the factor exp(— g^^^/2) accounts for the 
orbital character of Landau level wavefunctions. Since 
A{fio) is proportional to r~^, Eq. ( [T4|) must be solved 
self-consistentlv[15j giving rise to the following implicit 
equation for f = r/{e'^/£): 
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In Eq.([T7l) Aq = TA{na) and Vs = 2Trfns is the 'filling 
factor' of scatterers. Note that since the right hand side 
of Eq. ( I17p is a monotonically increasing function of Vg 
and a monotonically decreasing function of F, F must 
increase monotonically with i>s- 

Graphene QHF Phase Boundary — The Stoner criterion 
can be written in terms of F, Aq, and the dimensionless 
exchange integral 
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In the absence of screening (large F), X 
n = and X = (11/16) y^^ for n = 1. 
criterion, 

xAo/t > 1, 



= y^ for 
The Stoner 



(19) 



is satisfied for Vs > vl- Since v^ oc Us/B ex l/fiB, 
our Stoner criterion specifies the minimum values for the 
product of field and mobility illustrated in Fig. [TJ 

B [lOTesla] ^i [lO'^cm^V^^s^^] > 1/vI{vt). (20) 

The origin of the weaker tendency to ordered states in the 
four-fold degenerate n = 1 Landau level is the difference 
in form factor F{q). 

Discussion — Although the Stoner criterion can be ap- 
plied at all filling factors and provides a reasonable as- 
sessment of the crossover between interaction dominated 
and disorder dominated physics, we caution that the sim- 
ple quantum Hall ferromagnet states can occur only at 
integer values of the total filling factor. We expect the 
emergence of interaction-driven gaps at intermediate in- 
teger filling factors to be the first signal that sample 
quality is adequate to see interaction dominated physics. 
Judging by the relative size of charge gaps at integer 
and fractional filling factors, we expect that the first 
fractional filling factors will require mobilities approxi- 
mately five times larger than those required to realize 
quantum Hall ferromagnetism; the SU{A) nature of these 
Landau levels will open up a new frontier for the frac- 
tional quantum Hall effect that is likely to yield some 
surprises. We have so far neglected the Zeeman energy 
because it is much weaker than the competing disorder 
and interaction energy scales and will have little influ- 
ence on whether or not quantum Hall ferromagnetism 
occurs. (The Zeeman energy is ~ ImeV at ~ lOTesla 
compared to a n = Landau level interaction energy 
scale ~ lOOmeV depending on the degree of substrate di- 
electric screening.) When quantum Hall ferromagnetism 
does occur, however, the Zeeman energy will play a larger 
role. For v =^ ±1, in particular, the Zeeman energy will 
select ordered states that are spin-polarized, and break 
symmetry in the SU{2) valley space. The interaction 
terms in graphene should be weakly dependent on valley 
index, because interactions on the same graphene sublat- 
tice should be more strongly repulsive than interactions 
between sublattices at short distances, reducing the bro- 
ken symmetry to U{1). For this reason, we anticipate 
that the v — ±1 quantum Hall ferromagnet in graphene 
should have a Kosterlitz-Thouless phase transition at a 
low temperature. Finally we compare our result with 
the recent experiment [6| that has reported quantum Hall 
ferromagnetism in graphene. The mobility of the sample 
used in [6i] is /i = 5 X 10^[cm^/Fs]. Fig.[[T] indicates 
that for this mobility and v ~ ±1 the symmetry breaks 
at 17T in agreement with experiment [6[. The appear- 
ance of quantum Hall plateaus observed at v = ±4, in 



the middle of the four- fold degenerate n = ±1 Landau 
levels at around SOTesla is also in reasonable agreement 
with Fig.[[T] giving the critical field 40T. The infiuence 
of dielectric screening on our phase diagram, which we 
expect to be rather weak, and of screening due to vir- 
tual inter-Landau-level transitions, which we expect to 
be important at larger n, will be discussed in subsequent 
work. 

Acknowledgments-The authors acknowledge helpful in- 
teractions with Luis Brey, Andre Geim, Jason Hill, Zhi- 
gang Jiang, Philip Kim, Antonio Neto, Nuno Peres, Niko- 
lai Sinitsyn, and Jairo Sinova. This work has been sup- 
ported by the Welch Foundation and by the Department 
of Energy under grant DE-FG03-02ER45958. 



[1] 



[2] 



D. P. DiVincenzo and E. J. Mele, Phys. Rev. B 29, 1685 
(1984); C. L. Kane and E. J. Mele, Phys. Rev. Lett. 78, 
1932 (1997). 

G.W. Semenoflt, Phys. Rev. Lett. 53, 2449 (1984); F.D.M. 
Haldane, Phys. Rev. Lett. 61, 2015 (1988); C.L. Kane 
and E.J. Mele, Phys. Rev. Lett. 95, 146802 (2005); C.L. 
Kane and E.J. Mele, Phys. Rev. Lett. 95, 226801 (2005); 
N. M. R. Peres, F. Guinea, A. H. Castro Neto, Phys. Rev. 
B 73, 125411 (2006). V.P. Gusynin and S.G. Sharapov, 
Phys. Rev. Lett. 95, 146801 (2005); N.A. Sinitsyn et al, 
cond-mat/0602598 
[3] C. de C. Chamon et al, Phys. Rev. Lett. 77, 4194 (1996); 
D.H. Kim et al, Phys. Rev. Lett. 79, 2109 (1997); Y. 
Hatsugai and P. A. Lee, Phys. Rev. B 48, 4204 (1993); 

D. V. Khveshchenko, Phys. Rev. Lett. 87, 206401 (2001); 
J. Gonzalez et al, Phys. Rev. B 63, 134421 (2001); T. 
Stauber et al, Phys. Rev. B 71, 041406 (2005). 
Novoselov, K. S. et al. Science 306, 666 (2004); Y. Zhang 
et al, Phys. Rev. Lett. 94, 176803 (2005); C. Berger et 
al, J. Phys. Chem B 108, 19912 (2004). 
K.S. Novoselov et al. Nature 438, 197 (2005); Y.B. 
Zhang et al. Nature 438, 201 (2005). 
As this text was being written we learned that the antici- 
pated integer quantum Hall effects have been discovered; 
Y. Zhang et aL,|cond-mat/0602649 

Tomas Jungwirth and A.H. MacDonald, Phys. Rev. B 
63, 035305 (2001). 

D.P. Arovas et al, Phys. Rev. B 59, 13147 (1999); Z.F. 
Ezawa and K. Hasebe, Phys. Rev. B 65, 075311 (2002). 
A. A. Burkov, A.H. MacDonald, Phys. Rev. B 66, 115320 
(2002). 

T. Ando and Y. Uemura, J. Phys. Soc. Jpn. 37, 1044 
(1974); T. Ando et al. Rev. Mod. Phys. 54, 437 (1982); 
A.H. MacDonald et al, Phys. Rev. B 34, 2681 (1986). 

E. Stoner, Proc. R. Soc. London, Ser. A 165, 372 (1938). 
M. M. Fogler and B. L Shklovskii, Phys. Rev. B 52, 17366 
(1995); B.A. Plot et al, Phys. Rev. B 72, 245325 (2005). 
E. Fradkin, Phys. Rev. B 33, 3263 (1986); P. A. Lee, 
Phys. Rev. Lett. 71, 1887 (1993); K. Ziegler, Phys. Rev. 
Lett. 80, 3113 (1998); Y. Zheng and T. Ando, Phys. Rev. 
B 65, 245420 (2002). 

[14] S.V. Morozov, K.S. Novoselov, M.I. Katsnelson, F. 
Schedin, D. Jiang, A.K. Geim, |cond-mat/0603826, A. 



[4; 

[6; 
[7: 

[9; 

[10 



[11 

[12; 

[13 



F. Morpurgo, F. Guinea, cond-mat/0603789 , J. Gonza- [15] S. Das Sarma, Phys. Rev. B 23, 4592 (1981); Y. Mu- 
lez, F. Guinea, M. A. H. Vozmediano, Phys. Rev. B, 63, rayama and T. Ando, Phys. Rev. B 35, 2252 (1987). 

134421 (2001). 



